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Quantum fluctuations in an ultrafast rotating Bose gas at zero temperature are investigated. 
We calculate the condensate density perturbatively to show that no condensate is present in the 
thermodynamic limit. The excitation from Gaussian fluctuations around the mean fleld solution 
causes infrared divergences in loop diagrams, nevertheless, in calculating the atom number density, 
the correlation functions and the free energy, we find the sum of the divergences in the same loop 
order vanishes and obtain finite physical quantities. The long-range correlation is explored and the 
, algebraic decay exponent for the single-particle correlation function is obtained. The atom number 

(N . density distribution is obtained at the one-loop level, which illustrates the quantum fluctuation 

pH ' effects to melt the mean field vortex-lattice. By the non-perturbative Gaussian variational method, 

, we locate the spinodal point of the vortex-lattice state. 
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I. INTRODUCTION 

a 

The appearance of vortex excitations in response to rotation is a characteristic feature of superfiuid [ll-Q ■ Since the 
discovery of Bose-Einstein condensation (BEC) in atomic gases [1-11] , much work has been devoted to the properties of 
^ ' rotating gaseous condensates in traps, and these developments have been reviewed in 0, the rotation frequency 

^ , il increases, more and more vortices occur, from a single one to several ones, then they form an Abrikosov lattice, 
qh' i.e., a triangular array with a surface density n.v = ^-^{m is the mass of the condensed atoms) jol-fTsl. In the frame 
^ ; , co-rotating with harmonic traps, rotation effects are described by a centrifugal force, which effectively reduces the 
• transverse harmonic trapping, and a Coriolis force which has the same mathematical structure as the Lorentz force 
^ \ that an electron experiences in a uniform magnetic field. In the fast rotation regime, when the centrifugal force 
I i almost cancels the transverse confinement, the energy levels of the single particle part of the Hamiltonian organize 
' ^ ' into Landau levels with spacing 2h£L, and the condensates expand radically, leading to a very dilute atom number 
density, which ensures the mean interaction energy smaller than Landau level spacing, so that the cold atoms are 
confined in the lowest Landau level (LLL) single particle orbit. 

For rotating bosonic atoms in the LLL, the filling fraction, i.e., the ratio of the number of atoms to the number of 
vortices, is the parameter controlling the nature of the system. At high fillin g fra ctions, the condensate is in the mean 
field quantum- Hall regime and forms an ordered vortex lattice ground state |19l424| . As the filling fraction decreases, 
there is a zero temperature phase transition from a triangular vortex-lattice to strongly correlated vortex-liquid [25| 
and the melting point is located by various approaches j25l - [27| to be approximately at the filling fraction 6 '-^ 10. 
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[ — . The experiment done by V. Schwcikhard etc. |28j created an ordered vortex lattice in the mean-field quantum-Hall 
CsJ ■ regime, and provided evidence that the elastic shear strength of the vortex lattice drops substantially as the BEC 
enters the mean field quantum-Hall regime. 

In the ultrafast rotation limit when the transverse confinement is exactly canceled by the centrifugal force, the 
condensates expand to be a two dimensional configuration and atoms are frozen in the lowest energy level in the z 
direction (assuming a strong confinement in the z direction). For such a two dimensional system, J. Sinova etc. (26j 
find that the solution to the Gross-Pitaevskii (GP) equation is an Abrikosov triangular vortex-lattice and the integral 
• ^ • for the fraction of atoms outside the condensate diverges logarithmically with the system size, which implies that no 
' BEC occurs in the thermodynamic limit even at zero temperature. 

Fluctuation effects in the transverse plane in an ultrafast rotation limit are important and one has to go beyond 
mean field treatment p^ . In the present paper, we are devoted to quantum fluctuations in the ultrafast rotation limit 
and focus on the thermodynamic limit at zero temperature. 

Based on analytical calculations in the perturbative framework, we show that the condensate density is zero in the 
thermodynamic limit, namely, there is no BEC. The single-particle correlation function and the density fluctuation 
correlation function arc shown to fall off as an inverse power of the separation distance in the large distance limit, 
which indicates an algebraic long-range order and the algebraic decay exponent is obtained. The atom number density 
distribution is obtained at the one-loop level, which illustrates the quantum fluctuation effects to melt the mean field 
vortex-lattice. 

By loop expansion around the mean field solution, we calculate the free energy density up to two-loop. The mean 
field solution is an Abrikosov triangular vortex-lattice and the excitation from Gaussian fluctuations has a quadratic 
dispersion at small wave vectors [26| . We find that the quadratic dispersion causes infrared divergences in the two- loop 
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diagrams, nevertheless, the sum of the divergences vanishes and the two-loop contribution to the free energy density 
is finite. 

We also study the model by non-perturbative Gaussian variational method. The free energy density calculated by 
the perturbation theory and that by the Gaussian variational method coincide very well at large filling fractions. The 
vortex lattice solution exists only when the filling fraction v is greater than a certain value Vg-, which is found to be 
about 1.1. The point v = Vs vs, the so called spinodal point. Between the spinodal point and the melting point is the 
mcta-stable vortex-lattice state. 

The paper is organized as follows: 

In section [III we formulate the model. In section IIIIl we explore the long-range correlations. In section IIVI we 
calculate the free energy density perturbatively up to two-loop. In section |Vl we use the non-perturbative Gaussian 
variational method to study the model. In section IVll we give a summary and the conclusions. 



II. THE MODEL 



For a system of bosonic atoms in an axisymmetric harmonic trap (with trap frequencies wj^and ujz) rotating with 
angular velocity fie^ , the Hamiltonian in the rotating frame is 
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— is the strength of the hard core repulsive interactions, with Ug the s-wave scattering length. The 



where g 

centrifugal force effectively reduces the radial confinement and the Coriolis force is equivalent to the Lorentz force 
exerted on a particle with charge Q by a magnetic field B = ^^^e^,. In the present paper we are confined to the 
ultrafast rotation limit by setting D, = lo_i_, and assume the axial confinement is so strong that atoms are frozen in 
the lowest harmonic state in the z direction. Consequently, what we concern is essentially a two dimensional system 
of charged bosonic atoms experiencing an effective magnetic field in the z direction, described by the Hamiltonian 



N 



(pi — mflz X Tj)^ 



i=l 
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(2) 



The kinetic part of the Hamiltonian has equally spaced Landau levels (with spacing 2hfl) and the interaction part is 
a small perturbation, ng ^ 2hn (n is the mean number density of the atoms). Without thermal fiuctuations at zero 
temperature, all atoms are confined in the LLL. In the LLL subspace, the kinetic part of the Hamiltonian is quenched 
and represented by Hfl, hence we have the grand Hamiltonian in second quantized form (with the unit h = 1) 



{n - ti)¥{r)^{r) + |*^(r)^'^(r)*(r)*(r) 



H - fiN = J (fr 

and the grand-canonical partition function in the functional formalism 



(3) 



(4) 



where the action 4*] takes the form 



dT j Sy ^*{Y,T){dr- H + n)-^{v,T) + -g\^{T,T)\'^ 

in which (3 = -^^Kp and /i is the chemical potential. By variable rescaling : (^ — fi) 

^p^/3',r = --^=i=r', ^I* = ^/2mVl'^' , the partition function simplifies 



(5) 



Z'{f3',a^)^ j exp-^^ dr' j dV (r', r')(9.' - a^)vl''(r', r') + |*'(r', r')|^] 



(6) 



With all primes omitted, we have 



3 



Z{l3,a^,)^ Jv[^*^]exp-J^\r ^ d^r [^-*(r, T)(a. - a^)*(r, r) + |*(r, t)!^] . (7) 

Note that the kinetic part is absorbed in a shift of the effective chemical potential a^, and only the interaction part 
is relevant. Hereafter we will work with the rescaled variables. 

Using Landau gauge A = {—By,0) for the effective magnetic field B = ^^^e^, we have the LLL magnetic Bloch 
representation [2^, [s^l 



Vk(r) = 3» V exp 



— ;= — n Krd + ? — 71 — n + Zim— H 7iK„a + k^-x 

y/3\d 2 47r / l2^ ' d 2 ^ , 



(8) 



(p(r), i.e. <Pk(r) with k = 0, is a superposition of the lowest Landau levels of the kinetic part of the Hamiltonian and 
corresponds to an Abrikosov triangular vortex-lattice with lattice spacing d. In the rescaling length unit, d equals 

and the primitive vectors of the vortex-lattice are di = 0) , d2 = ^). describes the 

transverse oscillations of the vortex-lattice, and the primitive vectors of the reciprocal lattice are di = ^^(-^, — |), 
d2 = , /^(0, 1). The vortex cores, i.e. the zero point of the function ipk{r) with k = fcidi + A:2d2, are uniformly 



distributed at sites (ni — fc2)di + (712 + 5 + ki)d2 , where niand n2 are integers. Note that both the lattice cell and 
the Brillouin zone have an area of 2tt, and the vortex number density is n^, = 27r ^^^^ rescaling units. 

III. THE LONG-RANGE CORRELATIONS 

As is known j26| . for the ultrafast rotating two dimensional Bose gas, there is no BEC in the thermodynamic limit 
even at zero temperature. In spite of that, we will calculate U (1) invariant quantities, like the atom number density, 
the free energy and the correlation functions, in the perturbative framework, and show that the infrared divergences 
are canceled, similar to the method used in two dimensional non-linear a model [3l|, H^- In two dimensional 0{N) 



non-linear a model, David in [32| proved that, using the "wrong" spontaneously broken symmetry phase, any 0{N) 
invariant observable has an infrared finite weak coupling perturbative expansion. In this paper, though we will only 
show some U (!) invariant quantities are also free of infrared divergences at most to two loops in the perturbative 
framework, we believe that it is true to all orders similar to the two dimensional 0{N) non-linear a model. This 
method was extensively used to study the vortex lattice in type II superconductors, for example, in Ref. [33l.[33|. 

In this section. We calculate the condensate density perturbatively and show it is zero in the thermodynamic limit. 
By calculating the single-particle correlation function and the density fluctuation correlation function, we obtain the 
algebraic decay exponent. The atom number local density is also calculated, which shows that at large filling fractions, 
the number density retains the vortex-lattice configuration, while at small filling fractions, quantum fluctuations tend 
to smooth away the vortex-lattice. 

A. The atom number density 

In a usual fashion [sB - fss} . we separate the field as the condensate part and the fluctuation part 

«'(r,T) = V^^(r) +V'(r,T), (9) 

where Uc, the condensate number density, is a real number minimizing the free energy and the fluctuation part ■(/'(r, t) 
can be expanded as 

^{r, ^) = 7^ E E V'k™^k(r)e-3«'^e— . (10) 

^ P ]s.eBZ m 

In powers of V'km ^^"^ V'km, we divide the action S into four parts 
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So = l3A {-Qf^ric + PA-nf) , 

'52 = ^ [(-it^m - + 4ric/3k) V'kmV'km + "-c|7k| (V'kmV'-k-m + V'kmV'-k-m)] , 
P 

S'g = 2^^-^^ ^ (V'pi ^PlP2P30 +C.C.) , (11) 



P1,P2,P3 



= ]^ XI ^*pi^*p2^P3^P4Ppif 



Pi ,P2,P3,P4 

where p = (k, Wm) , V'p = ^km and J2p = Skesz • quadratic part can be diagonalized as 



5*2 = X {-ii^m + e(k)) flkn^akm , (12) 



p 

by the fohowing Bogohubov transformation: 



and inversely, 



where 



flkm = WkV'km + 1'kV'-k-m> (13) 
flkm = "kV'km +''-'kV'-k-m, (14) 



V'km = MkOkm - Ukfllk-m, (15) 
'/'km = Wkflkm - I'kfl-k-m, (16) 



2 V <(k) 



and 



2 V e(k) 



eo(k) = (-a^. + 4nc^k) , 



e(k) = V(-«A.+4«c/3k)'-4n2|7kP. (18) 



Note that the free energy density J- depends on two parameters, and Uc, nevertheless, ric is related to by the 



constraint '^-^("^■"a') _ q g^j^j hence ric is renormalized order by order. In the zero-loop order. 



J^c = -a^^^c + Pait-I, (19) 

and the condensate density equals 

= (20) 
2pA 

which is also the total number density in the zero- loop order, denoted as uq. Now we want to calculate the total 
number density and the condensate density to one-loop. The total number density n is given by 

1 J Sv {¥{v)^{v)) ^n^ + l^Y. (^kmV-k™) (21) 
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where ric is the condensate density and ^ J^p (V'kmV'km) is the density of the atoms outside the condensate. First 

we need to know the one-loop correction to the condensate density, denoted as ni^\ Up to one-loop order, the free 
energy density takes the form 

1 f cPk I 

TQ+i{af,,nc) = -af,nc + PAnl + — — ^ (-a^ + 47ic^k)^ - 4n2|7k|2 . (22) 



47r Jbz 27r 

Minimizing J-q+i with respect to ric leads to 



2/3^ 2nf3A Jbz 2n J(^-a,+AnAf-^n^c\l^? 



(23) 



from which we see that the one-loop correction to ric equals 



^(1)^ j rf2k 2/3k(2/3k-/3A)-|7k|\ (24) 



^"^^^^^ 2" v/(2/3k-/3A)^-|7k| 
The density of the atoms outside the condensate is given by 

where i5o(k) and £'(k) are defined as 

Eo{k) = 2/3k-/3A, (26) 



(25) 



Eik) = \/(2/3k-/3A)'-|7kP. (27) 

The total number density, rio+ii is equal to ric"' +ni^^ + (V'kmV'km)]^_;oop • The filling fraction, n/n-u, is given 

by 27r?io-|-i at the one-loop level. 

Obviously, ni^-* and -j^ J2p {^km''Pkm.) i^i^op both contain infrared divergences, but in the sum the divergences are 
canceled and ni^'^ + (V'km^km)i_;oop , which is the one-loop correction to the total number density, equals 

ni = / T^^(k) ^ -0-023. (28) 

^ttPa Jbz 27r 

The condensate density, ric = ^c"' +nc^\ is infrared divergent, ni'^^ = c(l — ahiL), where c and a are positive 

constants and L is the system size. If we can calculate ric to all loops, ric — c (l — a In L -I- i(a InL)^ -t- • • • ) ~ 
cexp(— alnL) = cL~". When we take L — >■ oo, ric — > 0. Therefore, in the thermodynamic limit, there will be no 
condensate, as is also shown by J. Sinova etc. [2g|. For a finite system, there is a finite infrared cutoff ~ the 
condensate density, ric = cL~°' , will be finite. 
Up to one-loop, the local density n{r) is given by 



(vl/t(r)vl/(r)) = (n(0) + n(i))^* (r)^(r) + ^ J^^ (r)^k(r) , (29) 

where n^'c''^ ip* {r)ip(r) is the mean field local density, denoted as no(r). Obviously, in Eq. ([^^ nc^^'and the last term 
both contain infrared divergences, nevertheless, by arranging them properly we see the divergences are canceled and 
the one-loop correction to the total local density is obtained as 

(1) , 1 f d^^Eoik) 



n,v*{v)^{v) + ^ j^^ (^k(r)^k(r) - ¥'*(r)^(r)) , (30) 



in which the second term is free of divergences. 
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FIG. 1: The mean field atom number density distribution forms a perfect triangular lattice, described by the function |iy3(r)p. 
Quantum fluctuations tend to smooth the mean field vortex-lattice. As the filling fraction lowers, quantum fiuctuations increase 
and the vortex-lattice becomes more smooth. 

Note that the mean field local density, no(r) — ni^'^ ip* {r)ip{r) , forms a vortex-lattice, and the one-loop correction, 
ni(r) , includes quantum fluctuations. In order to explore quantum fluctuation effects on the atom number density 



configuration, we plot the mean field density distribution, 

njl(r)+ni(r) 



"^l^*"^ = </j* (r)(^(r), and the total density distribution, 

no+ni ' different filling fractions. From Fig. [1] one finds that quantum fluctuations tend to smooth the vortex- 
lattice. At large filling fractions, the number density still retains the vortex-lattice configuration, while at small filling 
fractions, the vortex-lattice is smoothed away. 

B. The single-particle correlation function 

Up to one-loop, the single-particle correlation function (vl/t(i-j)vl>(r2)) is equal to 

(vE't(ri)vE'(r2)) = (n^ + n^) ^*(ri)v5(r2) + (V'^(ri)^(r2)) 

= (nf+AnW)^*(ri)^(r2) + A(^t(r,)^(r2)), (31) 

where A is used to keep trace of the loop order and should be set to 1 in the end. Similar to Eq. ([23), the above 
equation can be arranged as 



V3*(ri)(^(r2) + A/(ri,r2) 



.(°)+An(^)+AJ-/ 

inJsz 2n i?(k) 

+ Am) (p*(ri)(^(r2) + A/(ri,r2) , (32) 

where /(ri,r2) is given by^ /g^ [Vk(i"i)<^k(r2) - ip* {ri)ip{r 2)] and free of divergences. 

Now we define the relative coordinate r = ri — r2, and the center coordinate R = We take R = for 

simplicity and analyze the large |r| limit. Based on numerical calculations, we find 



kl 



lim /(ri,r2) ~ -0.18 In |r|<^* (ri)(p(r2) 



(33) 



Therefore 



lim (*1"(ri)*(r2)) ~ f Ji^^^ + Am) (ri)</5(r2) - 0.18A In |r|(^* (ri)(^(r2) 

|r|— foo \ / 

= (n(0) + Am) vp*(ri)(^(r2) (1 - aln|r|) 



(34) 



where a 



'^■^^^ and to one loop order a w '^'\a)' = ^7757 • Of course this result is not physical if we only include 



one-loop correction, as limr_i.oo (^^(ri))5'(r2)) ^ ^00 . We shall include all order in perturbation theory to get the 
physical result. We argue, in all order, similar to the calculation in Ref. [3^ . 



7 



lim (*t(ri)*(r2)) oc (ri)^(r2) 1 - a' In |r| + -(a' In |r|)2 + ... (35) 

r— ^oo \ I J 

= '/'*(ri)¥'(r2) Irp" 

To one loop, a' = a = similar calculations can be done in usual BKT (Berezinsky, Kosterlitz and Thouless) 

phase transition systems, and the algebraic decay exponent can be obtained correctly [39l. l40j. 

One will wonder in usual BKT phase transition systems, the phase transition is continuous. In the vortex lattice 
phase, from Ec^. ([551) . the correlation decays algebraically, and the rotational symmetry is broken as the factor 
Lp* (jci)(p{v2) in Eq. (|35p is not rotationally invariant. In the vortex liquid phase, the correlation decays exponentially, 
and the rotational symmetry is unbroken. Therefore, the phase transition is the spontaneous breaking of the rotational 
symmetry. As in usual solid to liquid phase transition, the phase transition is a first order melting transition. 

C. The density fluctuation correlation function 

In the following, we will calculate the density fluctuation correlation to one-loop, 

(a(ri)a(r2)) = ([n(ri) - {n{v,))] [n{v2) - (n(r2))]) 

= (n(ri)n(r2))-(n(ri))(n(r2)). (36) 

As shown in subsection IIII A[ there is no condensate and thus the system retains the U(\) symmetry, therefore, only 
the C/(l) gauge invariant terms in the contractions remain 

(fi(ri)n(r2)) - (n(ri)) (n(r2)) 
- (*t(ri)*(r2)) (*(ri)^'t(r2)) + {■^Hri)-^{Yi)-^Hv2)-^{v2)) ^ . (37) 

The term(^''''(ri)^'(r2)) is investigated in the last subsection. In the largejri — r2| limit, the first term in Eq. ([37|) 
falls off as an inverse power of the separation distance |ri — r2| , 

hm (vl/t(r0vl/(r2))(«'(ri)vl/t(r2))~|^(ri)h^(r2)|'|ri-r2r^". (38) 

|ri— r2|— i-oo 

The second term in Eq. (jST)) is the two-body connected Green's function and it is hard to calculate non-perturbatively. 
But we speculate that this term will not alter the asymptotic behavior of {5h{Yi)5n{T2)) in the largejri — limit, 
and hence the density fiuctuation correlation function also decays algebraically with the distance in the large distance 
limit. The result will indicate, near the Brag peak, the structure function S {Q + fc) where Q belongs to the reciprocal 
lattice and k is small will have a scaling Imik^Q S {Q + k) gc \k\~'^^'~^°'\ 

Similar calculation can be done for vortex lattice in type II superconductors and the density fluctuation correlation 
will be shown to have similar behavior. For the vortex lattice in type II superconductors near the Brag peak, the 
structure function S {Q + k) will have a scaling limfc_j.o S {Q + k) cx \k\~^ (details will be published elsewhere). 



IV. THE FREE ENERGY DENSITY CALCULATION BY LOOP EXPANSION 



In this section, we shall calculate the free energy density by loop expansion up to two-loop and show the cancelation 
of infrared divergences. 



A. Mean-fleld contribution 



In the saddle point approximation, 



SSl^^ SSl^^ (39) 
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with the LLL constraint, one obtains [sl. l4l| 



^o{r)^J^v{r), (40) 



2/3. 



A 



where (p{r) — (y5k=o(r) : /^A = ^ J^^n (fr\ip{r)\'^ . Obviously, the saddle point approximation is equivalent to the mean 
field GP equation. We have the mean field contribution to the free energy density 



^o(a^)=-T^, (41) 



and the mean field contribution to the number density 



no(a.) = -^?°^'^-^. (42) 



B. One- loop correction 



Following the loop expansion procedure presented in [42|, we set 



*(r,r) = J^^(r)+V(r,r), (43) 

in which y^^j(y5(r) is the mean field part and ?/'(r,T) is the higher order corrections, and then expand in 
powers of ^/'*(r, r) and "(/"(r, t). In the magnetic Bloch representation, ip{r, r) can be expanded as 

V'(r,T) = ^ ^ 5]^k™^k(r)e-5«-e"--- (44) 

where A is the area of the sample, Um = is the bosonic Matsubara frequency, and 6'k is defined in Eq. pS)) . In 
powers of V'km ^^'^ V'km: we divide the action into four parts 



4/3a y ' 

1 a,, 



5*0 = A/3 

^2 - E[(-*^'» + 1^ (2/5k - /3A)j ^■li^V'k™ + i^|7k| (V'LnCk-™ + V^k^V'-k-™)], 

" V^AS E (V'piV'p2^P3^PlP2P30+C.C.) , (45) 
\ PA P pi,p2,p3 

•^4 = -J-^ V'piV'p2'/'P3V'p4^PlP2P3P4) 
P1,P2,P3,P4 

where p = (k, w™) , ^/-p = Vkm , Ep = Ekesz Em ^^'^ 



/3k 


27r„ 


'cell 


7k 








'cell 




7k 
l7k| 


? 


p 

-'^PlP2P3P4 




-m2, 7713+7714 



(46) 



^^k,(r)^L(r)^k3(r)^k4(r)e^(^'^^+^''^-^>^3-«.4). 

cell ^TT 

In the one-loop approximation, we keep only the quadratic part of the action and diagonalizc it as 
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5*2 = XI + '^('*)) "kmOkm , (47) 

P 

where 



6(k)==^V'(2/3k-/3A)'-|7kP, (48) 
by the following Bogoliubov transformation: 

flkm = WkV'km + Uk-^-k-m ) (49) 

flkm = "kV'km + UkV'-k-m , (50) 

and inversely, 

'/'km = Wkflkm - ^'kalk-nu (51) 

'/'km = ""kflkm ^ Wkfl-k-m, (52) 



where Uk = f^^i^ + l) A'k = - 0' -^o(k) and i;(k) are defined in Eqs. ([26]) and By Taylor 



2 \^ £;(k) ^ -^y ''-'k — A/ 2 £;(k 

expanding /3k and |7k| , we find the excitation e(k) has a quadratic dispersion at small wave vectors, i.e. limk-i-o ^(k) ^ 
/c^, which is consistent with previous results (26l. I43ll43| . The one-loop contribution to the free energy density,J^i(a^), 
takes the form 



1 \ r ^ 

-^■^In J I?[a*a]exp-X(-iw„ + e(k))ak,„ak„ 

- y 



A 
weBZ 



i6(k) + iln(f-e-'3^x 



(53) 



By setting the area A to infinity and the temperature T to zero, we obtain 



Ma,)^£^{E{k))^, (54) 

where (•••)k = Ibz^^^ means average over the Brillouin zone. The one-loop correction to the number density, 
ni(ap), is equal to 

which is consistent with the result obtained in the last section as in Eq. 

C. Two- loop correction 

By the Bogoliubov transformation shown in Eqs. (jSip and (j52p . we switch to the field aj^^, Okm, and write the 
cubic and quartic part as 



^3 = ^Jj^^= Ki'^P2a;MpiP2P^ - K1P2P3) + ap,ap,ap,{llp^^^p^ - Up,p,p,)] + cc. (56) 



and 



Si — ^ Aj3 ^ ^ [Q'Pi'^P2'^P3'^P4^P1P2— P3— P4'^ki'L'k2''^k3'f'k4 



2ap^ap2Q'p3Q>p4 (-Ppi— p4P2P3'^ki'^k2'^k3''^k4 ~l" -Pp2P3Pl — P4 '^ki ''^k2 '^k3 'Uk4 )] + C.C.} 

X! a*^a*^ap3ap4[4Ppi_p,p3_p2UkiWk3''^k2'^k4 (57) 



Ap 

Pl:P2,P3,Pi 



~t~-fpip2P3P4 '^ki '^k2'^k3 '^k4 ~t- -fp3P4pip2 '^ki '-'k2 '^ks '^k4] 






FIG. 2: The two-loop Feynman diagrams. 
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where ap = Okm and 



HpipsM ^^Q (^'o-piPaPs'WkaWkaUki + -Po-psPaPi^kiWkat'k, + -P0-P3P2P1 "ki^ka ^kj ) , 



P1P2P3 



and 



A; 

a: 



P1P2P3 



P1P2P3 



(-Pp3P2-piOWk2Wk3Uki + JpiP3-p20UkiWk3'«k2 + ^pi P2 -paOUki ^ka "ka ) , 



^0piP3— P2^ki''^k3'^k2 ~(" -Pop2P3— Pi ^k2''-'k3Uki + -P0p3P2Pl ^ki Uk2''^k3 ) 
P2P3PlO'Wki'Wk3Uk2 + pip3P20^k2^k3^ki + -Ppip2P30^ki'L'k2''-'k3 ■ 



The two- loop contribution to the free energy density, J'2(a/x), takes the form 



1 



In 



JV[a*,a] exp-{S2 + S3 + Si) 



J 2?[a* , a] exp — ^2 



2 — loop 



(58) 



(59) 



(60) 



where (• • •) denotes the sum of ah the connected Feynman diagrams with Gp = Ve(k) ^ propagator. The 
two-loop Feynman diagrams are depicted in Fig. [2] The contribution from the diagram "cx)" equals 



'.A2 X! [^-^Pi- 



^ki^k2 "1' 4-Pp^_pjp2— P2 ^ki ^ki 'i^k2^k2 ~^ '^PpiP2PlP2 

1 



= ((l7k|iikt^k)k)' + ^ (/3ki-k2 K + O « + 0)k„ 

the contribution from the diagram "0" equals 



4iWk2 + "kiWk2)] GpiGp2 (61) 

(62) 



12a„ 



E In;. 



P1,P2,P3 



Op 

1? 



Gpi Gp2 Gp^ 



< 



nkik2(-ki-k2) nkik2(-ki-k2 



Eiki) + Eik2) + E{{ki +k2)) 



>ki,k2; 



(63) 



and the contribution from the diagram "OO "equals 

8a t, 



/3A/32A2 



2A2 E (^PiP2Pi ^piP2Pi) 



1 



{{[M< + ^^1) - l7k|«kt'k])J 



P3P2P3 
2 



a: 



P3P2P3 



(64) 



'ki,k2 



We have set the temperature T to zero and the area A to infinity in the end. The notation 

Ibz ^2^ Ibz ^2t^' means average over the Brillouin zone, and(ki + k2) represents the reduced wave vector in the 
Brillouin zone. Obviously, ^^-^2(0/^) = 0, and hence the two-loop correction to the atom number density is zero, 
77.2(0^) = 0. By Taylor expansion, one can confirm that each of the three contributions in Eqs. l|62|) . (|63p and ((64)) 
has infrared divergences. However, all the divergences are exactly canceled if they are summed up. By numerical 
integration, we find that 



Ma,.) = 

It is amazing that it precisely amounts to zero, and the reason is still under investigations. 



(65) 
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D. The filling fraction 



Up to two-loop, the free energy density has been obtained as 



and the atom number density equals 



-Fo+i+2(a^) = + (Eik))^ , (66) 



The fiUing fraction, ?t.o+i+2/'t-u, is given by 

TT 1 

PA ^PA 

Using the filhng fraction as a parameter, we have the free energy density 



{E{k))^. (68) 



^o+i+2(.) + am)) J' , (69) 

and numericaUy it is ^"0+1+2(1^) = — 0.029j^^ + 0.00062. To quantify the significance of quantum fluctuations, we 
calculate the ratio |— , given by j^^qqI^q , from which one can infer that at very large filling fractions, the quantum 
fluctuation effect is negligible and the mean field theory is adequate; at very small filling fractions, the quantum 
correction is even greater than the mean field contribution and the mean field vortex-lattice will be unstable due to 
drastic fluctuations. The quantum melting of the vortex-lattice has been intensively studied [25l - [27| . In the next 
section, we will locate the spinodal point of the vortex-lattice state, i.e. the terminal point of the meta-stable vortex 
lattice. 

V. GAUSSIAN VARIATIONAL CALCULATION 

In this section we are going to study the model by the Gaussian variational method. As in section Hill we set 

*(r,T) = V^^(r) V(r,T), (70) 
where Uc is a real number given by minimizing the free energy and ip(r, r) is expanded as 



where OJ^„j = 0_k-?ri , ^km ~ ^-k-jn- In powers of Okm and Akm, we divide the action S'[^'*4'] into four parts 

'S'c,'S'2,5'3,S'4. 



p 





( EO 






p) 




E^j 





S2 = ^yAO-p A^,)[ Ta]{ / ]■ (72) 



= -ap + 4nc/3k + 2nc|7k|, 
E^ = -a,, + 4n,/3k-2ne|7k|. (73) 



where Op = 0km , Ap = Akm , and 



Following the standard Gaussian variational procedure |45l - t47l | , we first define the Gaussian variational kernel 
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Gaussian variation 
loop expansion 



Gaussian variation 
loop expansion 



FIG. 3: The free energy density T and the fihing fraction v are shown. The blue square symbol represents the data obtained 
by Gaussian variational method, and the red circle symbol represents the result of the perturbation theory up to one-loop. 



where and ej^ are real variational parameters. The grand- canonical partition function can be written as 

I?[OA]e-«^^^(52-G + S'3 + 54r, (75) 



and the free energy density, J-"[nc, C?], is given by 



1 

JA 



Sc- In V[0 yl]e-« - ((-^2 - G + ^3 + ^4)") 



(76) 



where ()q denotes the sum of aU the connected Feynman diagrams with G as a propagator. Truncated to the first 
order, the free energy density J'[nc^ G] takes the form 



1 



Sc-\n j V[0 A]e-^ + {82- G + S3 + 5*4)0 



(77) 



Minimizing J^[nc,G] with respect to £^,£^ and Uc leads to the fohowing coupled equations 



£1. — -C/i, 



2tt 





l7k| 



1 7k' 
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FIG. 4; The relative difference, that is the ratio of the result of the Gaussian variational method minus that of the perturbation 
theory to the result of the perturbation theory, is shown. The blue square symbol represents the relative difference of the filling 
fraction, and the red circle symbol represents that of the free energy density. 

which can be solved numericaUy. Note that the effective chemical potential is the only physical parameter. 
Beginning with = 10 and slowly lowering it, we find that these equations cease admitting a solution when < 0.47, 
which indicates that the system can no longer support a vortex-lattice solution and hence it is the spinodal point. 
The atom number density can be obtained as the partial derivative of the free energy density with respect to the 
effective chemical potential and then the filling fraction is easily calculated. At the spinodal point, the filling fraction 
equals 1.1, smaller than that at the quantum melting point j25l - l27| . In Fig. [3] the free energy density and the filling 
fraction obtained by the Gaussian variational method are shown, together with those obtained by the perturbation 
theory and expressed in Eqs. (pS)) and To take a closer look at the very small differences of the results obtained 
by the two approaches, the relative difference, that is the ratio of the result of the Gaussian variational method minus 
that of the perturbation theory to the result of the perturbation theory, is plotted in Fig. |4l From the two figures, we 
see that the differences are very small and the higher the filling fraction, the smaller the relative difference. It reveals 
the fact that in the vortex-lattice state the perturbation theory and the Gaussian variational method coincide very 
well, and as the filling fraction gets higher, the two methods get closer and finally both are identical to the mean field 
theory in the large filling fraction limit. 

VI. SUMMARY 

We calculate the condensate density perturbatively to show that no condensate is present in the thermodynamic 
limit. By calculating the single-particle correlation function and the density fiuctuation correlation function, we 
obtain the algebraic decay exponent. We calculate the free energy density to two- loop and show the cancelation of 
the two-loop infrared divergences. The atom number density distribution to one-loop is obtained, which illustrates 
the quantum fluctuation effects to smooth away the mean field vortex-lattice. By the non-perturbative Gaussian 
variational method, we locate the spinodal point of the vortex-lattice, where the filling fraction, Vs, is numerically 
obtained to be about 1.1, lower than the quantum melting point obtained by various approaches [25l427l |. Between 
the spinodal point and the melting point is the meta-stable vortex-lattice state. From Fig. [1] we find that at the 
one-loop level, near the melting point where the filling fraction equals 6 ~ 10, the atom number density still retains 
the lattice configuration, while near the spinodal point where the filling fraction is about 1.1, the lattice is almost 
smoothed away. 

In order to determine the melting point accurately, we shall obtain the free energy density of the vortex liquid phase 
in the future. The study of the vortex liquid will be our focus in the future studies. 
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